Many-particle localization by constructed disorder: enabling quantum computing with 

perpetually coupled qubits 



^ ■ 
O ■ 

o ■ 
(N : 

^ . 

<: 

oo ■ 



M.I. Dykman\ F.M. Izrailev^, L.F. Santos\ and M. Shapiro^ 
^Department of Physics and Astronomy, and ^Department of Mathematics, 
Michigan State University, East Lansing, MI 48824 
^ Instituto de Fisica, Universidad Autonoma de Puebla, Puebla 72570, Mexico 

(Dated: February 2, 2008) 

We demonstrate that, in a many-particle system, particles can be strongly confined to their sites. 
The localization is obtained by constructing a sequence of on-site energies that efficiently suppresses 
resonant hopping. The time during which a many-particle state remains strongly localized in an 
infinite chain can exceed the reciprocal hopping frequency by > 10^ already for a moderate energy 
bandwidth. The results show viability of quantum computers with time-independent qubit coupling. 

PACS numbers: 03.67.Lx,72.15.Rn,75.10.Pq,73.23.-b 



I 

O 

o 



> 
o 

O 

o 

-I— > 



I 

o 
o 



X 



Disorder-induced localization has been one of the cen- 
tral problems of condensed-matter physics 0] . The prob- 
lem is particularly challenging for many-body systems, 
where only a limited number of results has been obtained 
0- Recently it has also attracted much interest in the 
context of quantum computing. In many proposed phys- 
ical implementations of a quantum computer (QC) the 
qubit-qubit interaction is not turned off j3, S S S 13 • 
Generally, the interaction leads to hopping of excitations 
between the qubits. Preventing hopping is a prerequisite 
for quantum computation. To control a QC and mea- 
sure its state excitations must remain localized between 
operations. Several approaches to quantum computation 
with interacting qubits were proposed recently |9l ll(lllll |. 

In this paper we study strong on-site many-particle lo- 
calization. It implies that each particle (or excitation) is 
nearly completely confined to one site. This is a stronger 
condition than just exponential decay of the wave func- 
tion, and it is this condition that must be met in a QC. 

Strong on-site localization does not arise in a disor- 
dered many- par ticle system with bounded random on- 
site energies [l^ ■ Indeed, consider a state where particles 
occupy sites. For short-range hopping it is directly 
coupled to ~ other iV-particle states. With proba- 
bility cx N one of them will be in resonance with the 
initial state, provided the on-site energies are uniformly 
distributed over a finite-width band. For large TV this 
leads to state hybridization over time ^ where J is 
the intersite hopping integral (we set h ~ 1). 

Here we show that many-particle localization can be 
obtained in an infinite chain by constructing a sequence 
of on-site energies. For the proposed narrow-band se- 
quence, all many-particle states remain confined for a 
time that largely exceeds J~^. We find that stationary 
many-particle states in moderately long chains are also 
strongly localized. 

A QC allows both studying and using strong local- 
ization. Here, on-site excitation energies are interlevel 
distances of the qubits. They can often be individually 
controlled, which makes it possible to construct an arbi- 



trary energy sequence. However, since the qubit tuning 
range is limited, so should be the energy bandwidth. A 
smaller bandwidth is also desirable for a higher speed of 
quantum gate operations. 

To localize one particle, the difference between excita- 
tion energies on neighboring sites should be much larger 
than J. In addition, even for nearest neighbor coupling, 
the energies of remote sites should also differ to prevent 
virtual transitions via intermediate sites. However, the 
further away the sites are, the smaller their energy differ- 
ence can be. We use this to obtain strong on-site single- 
particle confinement for a bounded energy bandwidth. 

For many-particle localization one has to suppress not 
only single-particle, but also combined resonances, where 
several interacting excitations make a transition simulta- 
neously. There is no known way to eliminate all such 
resonances. However, the problem can be approached 
from a different point of view by looking at a lifetime 
of a localized state. The effective many-particle hopping 
integral quickly falls off with the increasing number of 
involved excitations and intermediate nonresonant sites, 
which gives the effective "order" of a transition. To ob- 
tain a desired lifetime it is sufficient to eliminate reso- 
nances up to a certain order. As we show, this can be 
done for an infinite system. 

We will consider a one-dimensional chain of S — 1/2 
spins in a magnetic field. This model also describes 
an array of qubits. The excitation energy of a qubit 
is the Zeeman energy of a spin. The qubit-qubit in- 
teraction is the exchange spin coupling. For many pro- 
posed reafizations of QC's 0, 0, IE Hi it has a form 
i Y^' JI^^S^S^, where n, m are spin sites, ^ — x,y,z are 
spin projections, and J^,^ = for the effective mag- 
netic field in the z-direction. The ID spin system can 
be mapped, via Jordan- Wigner transformation, onto a 
system of fermions. For nearest neighbor coupling, the 
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fermion Hamiltonian is 

+JAV ajjaj,+ia„+ia„. (1) 

^ — 'n 

Here, ajj,a„ are the fermion creation and annihilation 
operators. Presence of a fermion on site n corresponds 
to the nth spin (qubit) being in the excited state. The 
on-site energies e„ in Eq. are the Zeeman energies 
counted off from the characteristic central energy, J = 
Jnn+i is the hopping integral, and JA = 
fermion interaction energy; we set J, A > 

Localization of stationary states can be conveniently 
characterized by the inverse participation ratio (IPR), 
which shows over how many sites the wave function 
spreads. For an A^-particle eigenstate IV'ata) (A enumer- 
ates the eigenstates) it is given by 



Jnn+i is the 
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where |$„i...„. 



, |0) is an on-site iV-particle 



wave function (quantum register). 

For fully localized stationary states — 1. For 
delocalized states In\ ^ 1 [for an L-site chain Ijyx < 
L\/N\{L — N)l]. Strong localization corresponds to /jva 
being close to 1 for all states. 

Localization requires that the on-site energies £„ be 
tuned away from each other. For nearest neighbor cou- 
pling a natural first step is to separate e„'s into two 
subbands, for even and odd n, with the inter-subband 
distance h that significantly exceeds J. Then we fur- 
ther split each subband into two subbands to detune next 
nearest neighbors. Here the splitting can be smaller, be- 
cause next-nearest-neighbor hopping occurs via a nonres- 
onant site, and the effective hopping integral is ~ J^/h. 
The procedure of band splitting is continued, with higher- 
order splitting being smaller and smaller. 

A simple sequence of e„ that implements the above 
idea has the form 



71 > 1 (3) 



([•J is the integer part). The energies © are illustrated 
in Fig. nja). Besides the scaling factor h, they are char- 
acterized by one dimensionless parameter a < 1. For 
small a, the two major subbands have width « ah and 
are separated by a gap of width « h. The splitting of 
higher-order subbands are proportional to higher powers 
of a. For a > 0.4 all subbands overlap and the subband 
structure disappears. 

One can see from Eq. © and Fig.QIa) that sites with 
close energies are indeed spatially separated. Analytical 
estimates of the energy difference can be obtained for 
small a. We have |e„+m — £n\ ~ ^ for odd m and ah 



for odd m/2. In general, the larger is m the higher may 
be the order in a of the leading term in |£ri+m — £n\- 

It is important for localization that the sequence 
(PI has no simple symmetry. It is neither self-similar 
nor quasi-periodic (wh ich is another example of "con- 
structed" disorder [13 )• For analytical estimates it is 
essential that the coefficients at any given power a' are 
repeated with period 2{q + 1) 

A convenient characteristic of the on-site energy se- 
quence is the amplitude of a particle transition from site 
n to site n + m. To the lowest order in J it has the form 



nm 
J/ [2(e„ - Sn+k)] 



(4) 



It can be shown using some results from number theory 
that Kn {m) decays with m nearly exponentially |l4l | . For 
small a and large |m| we have 



Kn{m) 



-v\m\ 



iJ/2h) 



|m| 



(5) 



The decrement v depends on rt, m. However, it is limited 
to a narrow region around = 1 with 0.89 < v < 1.19, 
cf. Fig. njb). For estimates one can use v — 1, i.e., set 
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FIG. 1: (color online). Single-particle localization for the on- 
site energy sequence 0. (a) The energies Sn/ft for a = 0.3. 
(b) The decrement v of the a-dependence of the transi- 
tion amplitude Knirn) Q for m = 1000 as function of n. 
The dashed lines show the analytical limits on v. (c) The 
mean single-particle inverse participation ratio (Ji) vs. a for 
hj J — 20 and for different chain lengths L. The vertical 
dashed line shows the analytical estimate for the threshold 
of strong localization. The inset shows the maximal IPR, 
-^1 max = maxA Ii\ , demonstrating strong localization. 

Equation Q describes the tail of the transition ampli- 
tude for J/2ha ^ 1. Strong single-particle localization 
occurs for a 3> ath, where the threshold value of a is 
ath ~ J/2h. The condition ath ^ a < 0.4 can be satis- 
fied for a moderately large ratio of the energy bandwidth 
h to the hopping integral J. 
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Strong single-particle localization for /i/ J = 20, as evi- 
denced by /lA being very close to 1, is seen from Fig.^c). 
The data are obtained by diagonalizing the Hamiltonian 
for open chains with different numbers of sites L. 

In the limit a ^ the stationary single-particle states 
are sinusoidal, which gives (/i) « L/3, cf. Fig.^c) ((•) 
means averaging over the eigenstates) . As a increases, 
the bands are split into more and more subbands, and 
(/i) decreases. It sharply drops to « 1 in a narrow region, 
which can be conditionally associated with a smeared 
transition to strong localization. The center of the tran- 
sition region gives ath- It appears to be independent of 
the chain length L. The estimate ath — J /'i.h is in good 
agreement with the numerical data for different h/ J . 

When a ^ ath, all states are localized. The wave 
function tails are small and limited mostly to nearest 
neighbors. At its minimum over a for given h/J, for 
all states Iw — 1 « J"^ see inset in Fig. ^c). The 
agreement with the above estimate becomes better with 
increasing hj J . For a > 0.4, when the bands of e„ start 
overlapping, the IPR increases with a. 

The difference in the localization problems for many- 
particle and single-particle systems stems from the in- 
teraction term cx JA in the Hamiltonian If A = 
(the interaction between the underlying spins or qubits 
is of the XF-type), the above single-particle results ap- 
ply to the many-particle system. For nonzero A, on 
the other hand, (i) the energy levels are shifted depend- 
ing on the occupation of neighboring sites, potentially 
leading to many-particle resonances, and (ii) there occur 
interaction- induced many-particle transitions. 

To analyze many-particle effects, it is convenient to 
change from ajj , a„ to new creation and annihilation op- 
erators bl^,bn that diagonalize the single-particle part of 
the Hamiltonian a„ = ^j^Unkbk- The interaction 
part of the Hamiltonian becomes 

H, = JAj2Vk,kMblA2^ksbk„ (6) 

where the sum runs over fci,2,3,4, and Vfcjfcafcafci = 

If all single-particle stationary states are strongly lo- 
calized, the off-diagonal matrix elements Unk are small. 
They are determined by the decay of the wave functions 
and fall off exponentially, Unk ~ K^''~"^ for \k — n\ ^ 1. 
At the same time, the diagonal matrix element is Unn ~ 
1. Therefore the major terms in the matrix Vkik2k3k4 are 
those with x — 0, where 

X- — min(|fci - p\ + \ k2 - p - l \ + \k3 - p - 1\ + \k4 ~p\). 
p 

They lead to an energy shift oc JA for the many-particle 
states with occupied neighboring sites. 

The terms 14ife2fe3fc4 with x > lead to two-particle in- 
tersite transitions (fca, ^4) — > (fci, k2), and FJA plays the 
role of a two-particle hopping integral. The parameter k 



gives the number of intermediate steps involved in a tran- 
sition. The steps are counted off from the configuration 
where the particles occupy neighboring sites |14lj . For 
large k and a ^ ath we have V oc {J/2ha'^)^ <^ 1, i.e., 
many-particle hopping integrals are small and rapidly de- 
crease with the number of involved virtual steps. 

Numerical results on the many-particle IPR are shown 
in Fig. 12 We have studied chains of length L = 10, 12, 
and 14 with L/2 excitations, which have the largest num- 
ber of states for given L {oc 2^ for large L). The results 
were similar, and we present the data for L = 12, in 
which case the total number of states is 924. 

For small a, the IPR is independent of a and is large 
because of the large number of resonating on-site states 
\^ni...ne)- It is reduced by the interaction cx JA that 
splits the energy spectrum into subbands depending on 
the number of occupied neighboring sites. On the whole, 
the IPR decreases with increasing a as long as a < 0.4. 
In the region 0.2 ^ a < 0.4 we have (le) ~ 1.01 except for 
narrow peaks. This indicates that away from the peaks 
all stationary states are strongly localized, as confirmed 
by the data on Jemax = max^ 
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FIG. 2: (color) Many-particle localization for a chain of length 
L — 12 with 6 excitations. The data refer to the first 12 sites 
of the chain the reduced bandwidth is h/J = 20. The 
purple, red, green, and blue curves give IPR for the coupling 
parameter A = 0, 0.3, 1, and 3, respectively. The inset shows 
the maximal Iq. Sharp isolated peaks for A 7^ result from 
the hybridization of many-particle on-site states that are in 
resonance for the corresponding a. The peaks for A = are 
due to the boundaries. 

A distinctive feature of the many-particle IPR as func- 
tion of a are multiple resonant peaks, a part of which is 
resolved in Fig. [21 They occur when two on-site states 
resonate, which gives Jemax ^ 2. The strongest peaks 
come from two-particle resonances. They happen when 
the two-particle energy difference 

Se = |£fci + ek2 - £k3 - £k4 1 (7) 

is close to Af JA with M = 0, 1, 2. 

As we increase a starting from a = 0, pronounced 
peaks of (Iq) appear first for Se « sah w J A with 
s — 1,2. They are due to hybridization of pairs on sites 
(n, n + 1) and (n, n + 3) for s = 1, and (n, n + 1) and 
(n — 1, n -|- 2) for s = 2, for example = 2-transitions). 
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For larger a, resonances occur when sd^h « MJA 
with n > 2. Such resonances require more intermediate 
steps, with K > 4. The widths of the IPR peaks are smaU 
and are in good agreement with simple estimates based 
on Eq. © In between the peaks /gmax = 1-02 for 

0-2 < a < 0.4 and h/J = 20. 

A special role is played by two-particle resonances 
where Se ^ J for all a < 0.4. They emerge already 
for >f = 2-transitions (n,n + 1) <-> (ri — l,rt + 2). Here, 



if n and n + 2 are prime numbers, Se 



h is ex- 



tremely small for large n. Strong resonance occurs for all 
n = 6fc — 1, in which case Se/h oc with ^ > 4. Such 6e 
is unusually small for k — 2. More many-particle reso- 
nances happen for higher >c. For different sections of the 
chain © we found that they could increase (/g) up to 
1.15 between the peaks, for hj J = 20, 0.2 < a < 0.4, and 
A = 1. These resonances can be eliminated by modify- 
ing the sequence (PJ: —> £„ + a'h/2 for n — 6k. For 
appropriate a' ~ 0.1, this modification brings (Jr) and 
-femax back to w 1.01 and « 1.02, respectively [lj,ll3- 

We now outline an alternative approach to the problem 
of strong localization, which is particularly relevant for 
quantum computing. Qubit states have finite coherence 
time, estimated as < 10^ — lO^J^^ for most models. It 
is sufficient to show that the states remain localized for 
a time that exceeds the coherence time. This will en- 
able both gate operations, that take time ~ J^^, and 
measurement, that often requires more time. 

Delocalization occurs through a transition to a resonat- 
ing on-site state. All resonant two-particle transitions 
up to a given number of steps >cq will be eliminated if, 
for ><• < xo, the on-site energy difference exceeds the 
maximum change of the interaction energy oc J A. For 
= 4 this requires JA/h < a^, a' /2, \a — ct' /2\. From 
Eq. ((BJ, the time needed for a transition with kq steps 
is > K^^'o/JA, it scales as h^'o . Then if ^ 4, the 
lifetime of all on-site states exceeds by a factor 10^ 
for h = 30,a = 0.3, a' « 0.2, and A < 1. 
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FIG. 3: (color online) The low-energy part of the two-particle 
energy differences Ssn/h O for all transitions with x < 5 in 
which one of the involved particles is on the nth site (n > 2). 
The data refer to a = 0.25. The left panel corresponds to the 
sequence Q. The right panel refers to the modified sequence 
with a' = 0.22 and shows the zero-energy gap. 

The transition time is further dramatically increased 
if >f = 4 resonances are eliminated. The low-energy gap 



in Se for the corresponding transitions is a^h when 
a' 3> a'^- For example, for a = 0.25 and a' = 0.22 this 
gap is 6e/h > 0.01 for all sites n > 4, cf. Fig. |3| This 
means that all states will remain localized on their sites 
for a time > K'^{JA)-^, if 2JA/h < 0.01. 

In terms of the operation of a quantum computer, the 
on-site energy sequence © is advantageous, since one 
radiation frequency can be used to resonantly excite dif- 
ferent qubits. Selective tuning to this frequency can be 
done without bringing neighboring qubits in resonance 
with each other. In our approach, localization does not 
require refocusing |^ , which is not always easy to imple- 
ment. We avoid delocalization due to indirect resonant 
n —> n + 2 transitions, which undermines the approach 
[lo| . Compared to Ref. |^, in our approach the interac- 
tion does not have to be ever turned off, and no multi- 
qubit encoding is necessary. In addition, our results are 
not limited to systems with nearest neighbor coupling. 

In conclusion, we have proposed a sequence of on-site 
energies that leads to strong localization of single- and 
many-particle stationary states of interacting spins or 
qubits. For an infinite chain, we eliminate resonances be- 
tween on-site states with the effective interstate hopping 
integral up to [J/h)^ . This leads to a long lifetime of lo- 
calized many-particle states already for a comparatively 
narrow energy bandwidth. The results show viability of 
scalable quantum computers where the interqubit inter- 
action is not turned off. 
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